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INTRODUCTION 


The  objective  of  the  project  is  to  test  whether  the  therapeutic  effectiveness  of 
permanent  implant  brachytherapy  for  prostate  cancer  can  be  improved  by  using  a 
combination  of  short  and  long  half  life  radionuclides  simultaneously.  Specific  aims  of 
the  proposed  project  are: 

1 .  To  test  theoretically  the  potential  of  a  mixture  of  radionuclides  in  permanent 
implants,  using  the  linear  quadratic  model,  as  a  function  of  Tpot,  potential  tumor  doubling 
time. 

2.  To  test  experimentally  the  validity  of  this  concept  by  in  vitro  irradiation  of 
BAl  112  sarcoma  cells  at  a  continuous  low  dose  rate  (CLDR)  with  (60  d  half  life), 
*®^Pd  (17  d  half  life)  and  a  50:50  mixture  of  and  *°^Pd  under  aerobic  conditions 
leading  to  exponential  growth  at  different  rates  (from  near  quiescence  to  full  exponential 
growth  at  a  maximal  rate,  with  a  doubling  time  of  approximately  14  hours). 

3.  To  measure  the  radiobiology  parameters  such  as  alpha,  beta,  half  life  of  repair  for 
the  BAl  1 12  sarcoma  cells  under  different  growth  conditions  and  develop  a  theoretical 
model  to  predict  expected  levels  of  cell  killing  using  *”^Pd  or  a  mixture  of  these 
isotopes. 

4.  To  use  immunohistochemical  techniques  to  measure,  in  solid  BAl  112  tumors  in 
vivo,  the  proportion  of  cells  in  S  phase,  the  proportion  proliferating  and  non-proliferating 
cells  and  the  tumor  doubling  time. 

5.  To  test  the  therapeutic  effectiveness  of  '®^Pd,  and  a  Pd/I  mixture  in  the 

BAl  1 12  in  vivo  tumor  system; 

6.  To  test  the  therapeutic  effectiveness  of  “^^Pd,  and  a  Pd/I  mixture  in  human 
prostate  carcinoma  xenografts  in  nude  mice,  using  a  slow  growing  and  a  fast  growing 
carcinoma. 

7.  To  evaluate  the  clinical  potential  and  feasibility  of  this  approach  in  the  treatment 
of  human  prostate  cancer. 

BODY  OF  THE  REPORT 

We  have  developed  a  theoretical  radiobiology  model  for  cell-killing  by 
continuous  low  dose  rate  irradiation  (CLDRI)  using  a  mixture  of  radionuclides. 
Theoretical  studies  were  performed  to  investigate  the  hypothesis  and  to  plan  in  vitro  and 
animal  studies.  Experiments  have  been  performed  using  BAl  1 12  tumor  cells  and  Chinese 
Hamster  cells  growing  in  vitro  and  BAl  112  cells  growing  in  vivo  as  solid  tumors  in 
WAG/rij  rats.  Radiobiology  parameters  for  these  cells  have  been  determined  and  used  in 
the  theoretical  radiobiology  model  to  improve  our  understanding  of  the  experimental 
observations. 
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The  linear-quadratic  model  of  cell-killing  by  CLDRI  Using  a  Mixture  of 
Radionuclides 


We  have  developed  a  theoretical  model  for  CLDRI  using  a  mixture  of 
radionuclides  with  different  half  lives.  This  model  is  described  in  the  attached  manuscript 
entitled  “Biologically  effective  dose  (BED)  for  interstitial  seed  implants  containing  a 
mixture  of  radio-nuclides  with  different  half  lives”  by  Zhe  Chen  and  Ravinder  Nath. 
Briefly,  the  purpose  of  this  project  was  to  develop  a  tool  for  evaluating  interstitial  seed 
implants  that  contain  a  mixture  of  radionuclides  with  different  half-lives  and  to  examine 
the  clinical  implications  of  prescribing  to  an  isodose  surface  for  such  an  implant.  Using  a 
generalized  equation  for  the  biological  effective  dose  (BED)’"^,  the  effects  of  cell 
proliferation  and  sub-lethal  damage  repair  were  examined  systematically  for  implants 

containing  a  mixture 
of  radionuclides 
(Figure  1).  The 
results  were 
contrasted  with  those 
for  implants  using  a 
single  type  of 

50  radionuclide.  A 

40  BED  (Gy)  clinical  permanent 

seed  implant  that 
contained  a  mixture 
of  and  ’°^Pd  seeds 
was  used  to  examine 
the  clinical 
implications  of  the 
isodose  prescription 
for  such  implants.  An 
equation  of  BED  for 
implants  containing 
any  number  of 
radionuclide  types 

was  obtained.  For  implants  containing  a  mixture  of  radionuclides  with  different  half-lives 
such  as  and  *®^Pd,  the  dose  as  well  as  its  temporal  delivery  pattern  to  a  point  is 
dependent  on  the  relative  dose  contributions  from  different  types  of  radionuclide.  It  can 
vary  from  point  to  point  throughout  the  implant  volume.  Therefore  the  quantitative 
effects  of  cell  proliferation  and  sub-lethal  damage  repair  are  spatially  dependent  in  such 
an  implant.  For  implants  containing  a  mixture  of  ’^^I  and  *”^Pd  seeds,  the  prescription  to 
an  isodose  surface  becomes  non-unique.  If  the  prescription  dose  was  based  on  existing 
clinical  experience  of  using  *^^I  seeds  alone,  mixing  *  ^Pd  seeds  with  '^^I  seeds  would 
decrease  the  cell  survival  in  such  implant.  On  the  other  hand,  if  the  prescription  dose  was 
based  on  existing  clinical  experience  of  using  *®^Pd  seeds  alone,  mixing  I  seeds  with 
*®^Pd  seeds  in  a  same  implant  would  create  radiobiologically  “cold”  spots  (i.e.  an  increase 
in  cell  survival  from  the  clinical  expectation)  at  locations  where  a  major  portion  of 
prescription  dose  is  contributed  by  the  *^^I  seeds.  For  fast-growing  tumors,  these  "cold" 
spots  can  become  significant.  From  this  theoretical  investigation,  we  conclude  that  when 


D(’^®l)/TD 


Figure  1 
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cell  proliferation  and  sub-lethal  damage  repair  are  present  during  dose  delivery,  total  dose 
alone  is  no  longer  sufficient  for  a  complete  characterization  of  an  interstitial  seed  implant. 
In  order  to  avoid  radiobiological  “cold”  spots  when  radionuclides  of  different  half-lives 
are  mixed  in  a  permanent  implant,  the  dose  prescription  should  be  based  on  the  clinical 
experience  of  using  the  longer  half-life  radionuclide.  Biologically  effective  dose  provides 
a  tool  to  start  examining  the  radiobiological  effects  of  mixing  different  type  of 
radionuclides  in  the  same  implant. 

A  manuscript  has  been  submitted  for  publication  in  the  International  Journal  of  Radiation 
Oncology. 


In  vitro  CLDRI  studies 


BAl  112  tumors  were  grown  between  the  ears  of  the  a  14  week  old  male  WAG/rij 
rats  by  interdermal  inoculation  from  a  single  cell  suspension  of  BAl  112  cells  obtained 
from  a  21  day  BAl  112  tumor  growing  on  the  head  of  a  previously  inoculated  rat.  A 
tumor  cell  suspension  was  made  from  the  BAl  1 12  tumor  and  between  1.5x10^  and 
5.0x10^  cells  were  plated  into  petri  dishes.  These  cells  were  allowed  to  settle  and  reach 
logarithmic  growth,  (48-72  hrs),  before 
they  were  used  in  a  continuous  low 
dose  rate  experiment. 

Monolayers  of  rat 
rhabdomyosarcoma  cells  (BAl  112) 
were  irradiated  in  vitro  by  '°^Pd  sources 
in  a  polystyrene  phantom.  Colony 
formation  ability  of  irradiated  cells 
under  aerobic  conditions  was  measured 
for  graded  doses,  at  a  dose  rate  of  6 
cGy/hr.  Dose  to  the  cell  monolayers 
was  determined  using  FeS04  Fricke 
dosimetry,  with  a  calculated  correction 
for  interface  effects  due  to  photoelectric 
effect  in  the  tissue  culture  dishes.  The 
sources  (up  to  60  in  one  experiment) 
were  arranged  in  concentric  circles  in 
such  a  way  as  to  provide  a  dose 
uniformity  of  better  than  ±  5%  across 
the  dishes.  Some  of  the  results  are 
shown  in  Figure  2.  Comparison  of  the 


Figure  2. 


surviving  fraction  as  a  function  of  dose  calculated  for  the  BAl  112  cells  to  that  measured 
by  using  CLDRI  *°^Pd  irradiation  is  also  shown  in  Figure  2.  The  lines  through  the  data 
points  represent  the  calculated  survival  curves  the  symbols  represent  the  measured  data 
The  parameters  of  a,  p,  repair  half-time,  and  tumor  doubling  time  were  determined 
directly  from  the  measurements  performed  on  the  BAl  112  cells,  as  described  later  in  the 
report. 
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In  vitro  studies  at  an  acute  dose  rate  using  simulated  x-rav  beams 


To  study  the  radiobiological  characteristics  of  the  cells  under  acute  exposure 
condition,  simulated  x-ray  beams  with  average  energies  equivalent  to  that  emitted  by 
(27.2  -  35.49  keV  with  an  average  of  27.4)  and  “’^Pd  (20  -  22.7  keV  with  an  average  of 
20.5  keV)  were  established  on  a  new  orthovoltage  unit.  The  simulated  beams  not  only 
have  the  average  energies  similar  to  that  given  by  the  radioactive  isotopes  but  also  have  a 
narrow  photon  energy  spectrum.  The  narrow  photon  energy  spectrum  was  achieved  by 
optimizing  the  tube  voltage  (which  determines  the  upper  limit  of  the  produced  photon 
energy)  and  the  added  filtration  (which  filters  out  the  low-energy  Bremsstrahlung 
photons),  following  the  work  of  Muench  et  al^. 

Aluminum  filters  fi-om  Pantak  was  used  to  construct  a  customized  filter  for  the 
DXT  300  unit  that  has  a  desired  filtration  thickness.  A  set  of  aluminum  filters  (with 
thickness  of  0.1  to  1.0  mm)  fi'om  Nuclear  Associates  (AL  Filter  Set  07-430)  were  used  to 
determining  the  half- value-layer  (HVL)  of  a  simulated  beam  using  a  customized 
filtration.  The  thickness  of  the  aluminum  sheets  were  measured  by  using  a  Mitutogo 
micrometer  (Serial  #  2032360)  with  accuracy  of  0.001  mm).  The  aluminum  HVL  for  a 
given  beam  was  determined  by  in-air  ionization  chamber  measurement  under  the  narrow 
beam  geometry.  An  air-equivalent  Spokas  chamber  (Exradin,  Model  No:  A1  (0.5  ml,  AE 
plastic))  was  used  to  measure  the  ionization  at  a  fixed  source  to  chamber  distance  (SCD) 
of  50  cm  in  air.  The  ionization  charge  was  measured  by  a  Keithley  electrometer  (model 
35614E  SN  43075)  with  -300  V  bias  potential.  Due  the  energy  dependence  of  the 
chamber  at  the  low  energies,  the  measured  ionization  were  converted  to  corresponding 
exposure  and  the  HVL  is  then  determined  from  the  relative  exposure  as  a  function  of 
aluminum  filter  thickness.  A  narrow  circular  beam,  with  a  diameter  of  6  cm  at  SCD  of  50 
cm,  was  generated  by  using  a  homemade  lead  collimator  mounted  to  DXT  300’s 
accessory  mount.  The  aluminum  sheets  added  to  the  beam  for  the  HVL  measurement 
were  taped  to  the  bottom  of  the  lead  collimator. 

The  HVL  as  a  fimction  of  mono-energetic  photon  beam  energy  for  aluminum  is 
taken  from  Johns  and  Cunningham^.  The  expected  HVL  for  a  simulated  (*®^Pd)  beam 
with  average  energy  of  27.4  keV  (20.5  keV)  is  1.84  mm  (0.82  mm)  aluminum.  With  the 
expected  HVL  in  mind,  the  tube  kV,  mA  and  the  thickness  of  the  added  filtration  were 
optimized  for  a  simulated  and  a  simulated  *”^Pd  x-ray  beam.  The  optimum  setting 
determined  for  the  DXT  300  unit  is  summarized  in  the  following  table  I. 


Table  I.  Radiation  characteristics  of  simulated  x-ray  beams 


Beam 

kV 

mA 

Added  Filter 
(mm  AL) 

Beam  HVL 
(mm  AL) 

Energy 

Homogeneity 

(%) 

Equivalent 

Energy 

(keV) 

<E>  from 
isotope 
(keV) 

1-125 

Equivalent 

43 

20 

3.545 

1.851 

86.9 

27.45 

27.4 

Pd-103 

Equivalent 

29 

25 

1.826 

0.82 

88.6 

20.5 

20.5 

Results  of  the  cell  survival  experiments  at  an  acute  dose  rate  using  these  x-ray 
beams  are  shown  in  Figure  3  and  4.  The  parameters  of  a,  p,  repair  half-time,  and  tumor 
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doubling  time  were  determined  directly  from  the  measurements  performed  on  the 
BAl  1 12  cells  (details  are  provided  separately). 


> 

'E 
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Figure  3.  Comparison  of  the  surviving 
fraction  as  a  function  of  dose  calculated  for 
the  BAl  112  cells  to  that  measured  using  the 
simulated  x-ray  beams. 


Figure  4.  Comparison  of  the  surviving 
fraction  as  a  function  of  dose  calculated  for 
the  BAl  112  cells  to  that  measured  using  the 
simulated  '“Pd  x-ray  beams. 


Procedure  for  in  vitro  studies  for  quiescent  BA1112  cells 

A  procedure  for  in  vitro  studies  for  quiescent  BAl  112  cells  has  been  developed.  For 
these  studies,  animals  will  be  implanted  with  transplanted  BAl  112  rhabdomyosarcomas 
by  inoculation,  into  a  subcutaneous  site  on  the  heads.  The  tumors  will  be  allowed  to 
grow  for  3  weeks,  to  an  experimental  volume  of  approximately  100-200mm^.  The 
animals  chosen  for  the  quiescent  cell  experiments  will  be  euthanized  by  anesthetic 
overdose  and  the  tumor  cells  will  then  be  removed  using  aseptie  techniques.  A  single¬ 
cell  suspension  of  tumor  cells  will  be  suspended,  counted,  and  assayed  for  viability  using 
the  same  colony  formation  assay  used  for  cells  in  cultures.  1.5x10^  cells  will  be  plated 
into  a  flask  with  13  ml  of  DMEM  for  cell  growth.  These  cells  will  then  be  passed  twice  a 
week  for  approximately  4-8  passages.  The  cells  are  transplanted  for  2-4  weeks  to  assure 
a  homogeneity  population  of  BAl  1 12  cells.  After  8-10  passages  of  the  BAl  1 12  cells,  in 
vitro,  a  new  in  vivo  tumor  is  used  to  produce  a  new  primary  single-cell  passage,  as 
described  above. 

When  the  first  line  has  reached  passage  5  the  cells  are  eounted  and  plated  into  5  petri 
dishes  for  a  growth  curve  study.  1.0  xlO^  cells/dish  are  planted  on  Day  0.  The  cells  are 
allowed  to  grow  for  24  hours  and  cell  counts  are  taken  on  days  3-7.  Our  first  cell  count 
reading  is  taken  24  hours  after  initial  seeding.  This  is  to  assure  lag  phase  has  occurred 
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and  the  cells  are  in  exponential  growth.  Within  this  one  week  exponential  cell  growth 
occurs,  and  a  population  of  3.6x10^  cells/dish  accumulates  by  Day  6.  Plateau  phase  is 
reached  by  Day  7  (3.1x10®  cells/dish).  An  in  vitro  growth  curve  is  plotted  to  assure  the 
primary  tumor  cells  have  adjusted  to  their  in  vitro  environment.  Their  doubling  times  are 
approximately  24  hours. 

In  our  initial  experiment,  the  exact  cell  number  per  dish  was  repeated  as  in  the  previous 
experiment  with  the  addition  of  5  dishes  (media-change  dishes).  The  media  in  these 
dishes  are  changed  every  24  hours,  to  accelerate  cell  division  and  to  reach  quiescent 
growth  more  quickly.  Upon  completion  of  this  experiment,  the  control  dishes  (media  not 
changed),  showed  similar  results  as  the  control  experimental  dishes,  but  the  experimental 
dishes  (media  replaced)  had  a  dramatically  higher  cell  counts  and  did  not  seem  to  reach  a 
plateau  phase.  The  experimental  dishes  reached  1.1x10^  cells/dish  and  Day  7,  and  were 
not  growing  exponentially  as  in  Days  3-5,  but  were  in  a  quiescent  phase,  still  alive  and 
not  dividing. 

In  our  third  experiment,  cells  were  plated  similarly  as  the  previous  experiment  but  the 
initial  seeding  of  the  dishes  was  increased  to  10®  cells/dish.  Our  results  showed  the 
control  dishes  seeded  with  10®  cells/dish  and  no  media  changes,  throughout,  reached 
plateau  by  Day  5  and  the  quiescent  phase  was  probably  attained  by  late  Day  3-4.  The 
experimental  dishes,  seeded  with  10®  cells/dish  also,  and  the  media  replaced  every  24 
hours,  showed  early  plateau  phase  on  Day  7.  Therefore,  the  quiescent  stage  was 
somewhere  between  Days  5-7. 

Both  protocols  will  be  useful  in  staging  of  cells  for  early  or  late  week  quiescent  cell 
survival  curve  experiments  on  the  Pantak  250  x-rays,  simulating  1-125  and  Pd- 103 
energy’s. 

Radiological  Parameters  for  CCL-16  cells 

In  order  to  measure  the  radiological  parameters 
of  the  cells  used  in  these  studies,  the  acute 
exposure  survival  curves  were  measured  using 
250  kV  x-rays.  Split  dose  experiments  were 
conducted  to  measure  the  half  time  of  sublethal 
damage  repair.  The  solid  line  in  Figure  5 
represents  the  measured  survival  curve  for 
CCL-16  cell  line  for  the  250  kV  x-rays.  Solid 
symbols  represent  acute  exposure  while  the 
open  symbols  represent  the  survival  when  the 
single  exposure  dose  is  split  into  two  equal 
doses  separated  by  2  hour  interval.  The  cells 
are  kept  at  22®C  during  the  waiting  interval  to 
prevent  proliferation.  The  a  and  P  parameters 
were  estimated  by  fitting  the  acute  exposure 
curve  to  a  LQ  equation.  The  a  and  p  for  the 
250  kV  acute  exposure  were  found  to  be  0.42 


Figure  5 
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Gy*  and 0.023  Gyl 

With  p  determined  from  the 
acute  exposure  survival  curve,  the  _ 
sub-lethal  damage  repair  time  constant  | 
can  be  determined  from  the  ratio  of  s 

the  survival  measured  from  the  split  5 

dose  irradiation  to  that  of  single 
irradiation.  The  logarithm  of  this  ratio 
is  a  quadratic  function  of  total  dose 
(Figure  6).  The  repair  half  time,  using 

the  incomplete  repair  model,  is  found  to  be 
1.1  hours. 
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Acute  survival  curves  for  BAl  112  cell  line  for  acute  exposures  using  the  simulated  I 

and  *°^Pd  x-ray  beams  were 
measured.  Results  are  shown  in 
Figure  7.  The  lines  represent  fit  to  LQ 
model  using  a  fixed  p.  The  parameters 
of  a  and  P  for  the  BAl  112  cell  line 
for  and  ’®^Pd  photons  were 
determined  from  the  survival  curves 
obtained  under  acute  exposure 
condition  using  simulated  and 
*®^Pd  x-ray  beam  on  an  orthovoltage 
unit.  The  acute  survival  curves  for  the 


simulated  and  ‘”^Pd  x-rays  were 
first  fitted  individually  to  LQ  equation 
for  a  and  p.  The  average  of  the  p 
from  these  fits  is  then  used  as  a  fixed 
parameter  in  the  fit  of  survival  curves 
to  determine  the  a  for  each  simulated 
x-rays.  The  fitted  a  and  p  is  given  in 
table  II. 
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Split  dose  experiments  were 
performed  with  250  kV  x-rays  for 
BAl  1 12  cell  line.  Results  are  shown  in  Figure  8.  A  total  dose  of  10  Gy  was  given.  The 
survival  fraction  was  measured  as  a  function  of  time  interval  between  two  5  Gy 
irradiation.  The  cells  were  kept  at  22”C  during  the  waiting  interval  to  prevent 
proliferation.  According  to  incomplete  repair  model  we  have. 


Figure  7 
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ln[^(0/5(0)] 

lii[X(oo)/S(0)] 

where  S(t)  is  the  surviving 
fraction  with  split  time 
interval  t  and  p  is  the  sub- 
lethal  damage  repair  time 
constant.  By  fitting  the 
measured  data  to  the  above 
equation,  the  repair  half-time 
was  found  to  be  19  minutes. 
The  following  table  II 
summarizes  the  parameters 
determined  for  the  BAl  112 
cell  lines. 

Tumor  growth  was 
measured  for  the  BAI 112  cell 


Split  Interval  (hours) 


Figure  8 


implanted  in  rat  (Figure  9).  The 

tumor  growth  curve  is  fitted  to  an  exponential  function  of  time.  The  tumor  size  doubling 
time  is  found  to  be  2.7  days. 
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Table  II.  Radiological  parameters  for  BA1112 


Parameters  for  BA1112 

Mean  photon  energy 

a 

(Gy') 

(Gy-') 

a/p 

(Gy) 

Tumor  Doubling 
Time 
(day) 

Sub-lethal  Damage 
Repair  Half-time 
(hr) 

equivalent 

0.35 

0.035 

10.1 

2.7 

0.314 

equivalent 

0.28 

0.035 

7.9 

2.7 

0.314 

In  vivo  tumor  growth  studies 

In  order  to  produce  a  consistent  dose  distribution  to  irradiate  the  tumors 
transplanted  to  different  animals  and  to  minimize  radiation  exposure  to  personnel 
handling  the  radioactive  seeds,  an  afterloading  seed  applicator  was  designed  and 
fabricated.  The  applicator  was  made  of  polystyrene  with  loading  ports  for  nine  seeds.  The 
central  portion  of  the  applicator  was  open  and  has  a  dimension  large  enough  for  tumor  to 
grow.  Equal  source  strength  was  assigned  to  all  nine  seeds  in  order  to  minimize  the 
possible  confusion  of  handling  variable  source  strengths.  The  seeding  configuration  was 
optimized  to  produce  an,  as  uniform  as  possible,  dose  distribution  to  the  central  portion  of 
the  applicator  and  to  be  usable  for  both^^^I  and  '®^Pd  seeds. 

Twelve  applicators  were  built  to  conduct  a  tumor  cure  experiment.  A  light¬ 
weight  metallic  helmet  measuring  2.2  cm  (cranial-caudal  length)  x  2.15  cm  (side-side 
width)  X  2.2  cm  (height)  was  sutured  to  the  rat’s  head  by  four  stitches  through  the 
cartilage  of  the  ears  and  two  more  stitches  just  behind  the  head  at  the  neck.  A  seventh 
suture  placed  under  the  tumor  will  be  tied  to  the  central  bar  that  across  the  top  of  the 
helmet,  thus  ensuring  the  tumor  is  pulled  up  into  the  center  of  the  treatment  volume. 
Treatment  volumes  were  2.3  cm^  for‘^^1  or  ‘®^Pd.  The  lightweight  metallic  helmet  was 
afterloaded  with  the  seed  applicator.  Calculated  dose  rates  from  each  applicator  are 
shown  in  tables  III  and  IV. 
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Table  III.  Relative  Initial  Dose  Rate  for  Draximage  LS-1  Seeds 


Dose  Points 

Relative  Dose  (%) 

Center 

3  mm  (X) 

4  mm  (X) 

5  mm  (X) 

3  mm  (D)  4  mm  (D) 

5  mm  (D) 

9  mm  up 

8  mm  up 

126.7 

115.1 

105.4 

116.9 

109.0 

112.7 

113.6 

119.1 

7  mm  up 

105.6 

109.4 

119.4 

117.2 

112.3 

117.4 

125.2 

6  mm  up 

106.0 

113.4 

121.0 

122.6 

112.7 

121.5 

133.9 

5  mm  up 

104.6 

114.7 

122.9 

133.0 

115.8 

127.5 

142.6 

4  mm  up 

103.1 

115.5 

131.6 

143.2 

120.3 

130.6 

145.5 

3  mm  up 

104.4 

119.2 

135.6 

148.9 

122.3 

132.9 

143.8 

2  mm  up 

103.4 

120.7 

135.5 

149.7 

121.2 

133.5 

140.5 

1  mm  up 

100.0 

119.0 

134.0 

147.0 

117.2 

130.0 

137.3 

Base 

94.8 

111.8 

125.9 

137.6 

111.2 

121.2 

127.6 

1  mm  down 

87.2 

99.7 

112.0 

121.1 

101.7 

110.7 

114.6 

2  mm  down 

78.1 

84.8 

94.4 

100.6 

89.3 

96.6 

100.1 

Table  IV.  Relative  Initial  Dose  Rate  for  *”Pd  Theragenic  M200  Seeds 


Dose  Points 

Relative  Dose  (%) 

Center 

3  mm  (X) 

4  mm  (X) 

5  mm  (X) 

3  mm  (D)  4  mm  (D) 

5  mm  (D) 

9  mm  up 

8  mm  up 

120.2 

114.9 

110.0 

105.9 

101.2 

115.3 

115.1 

113.6 

7  mm  up 

112.0 

110.1 

109.8 

107.5 

114.2 

118.1 

120.3 

6  mm  up 

107.3 

112.9 

116.1 

116.2 

115.9 

122.0 

126.4 

5  mm  up 

105.3 

116.7 

123.7 

125.8 

118.1 

126.4 

131.7 

4  mm  up 

106.5 

120.4 

129.0 

132.3 

119.8 

129.9 

134.6 

3  mm  up 

107.1 

122.8 

130.1 

135.6 

122.6 

131.6 

135.8 

2  mm  up 

104.4 

122.3 

130.3 

136.5 

122.4 

128.9 

131.8 

1  mm  up 

100.0 

119.2 

126.1 

132.7 

119.1 

124.4 

123.2 

Base 

94.9 

112.2 

118.9 

124.2 

111.4 

115.9 

113.5 

1  mm  down 

88.1 

100.9 

107.3 

111.1 

99.7 

103.8 

101.5 

2  mm  down 

78.9 

88.3 

93.1 

95.1 

84.7 

88.0 

86.7 

The  actual  dose  delivered  for  a  given  batch  of  seeds  were  verified  by  TLD  point 
dose  measurements  for  each  applicator  after  it  is  removed  from  a  rat  at  the  completion  of 
experiment.  A  jig,  which  models  the  fully-grown  tumor  on  the  head  of  a  rat,  is  made  of 
polystyrene.  Two  1x1x1  mm^  micro-TLD  cubes  can  be  placed  in  the  jig  so  that  the  TLD 
cube  is  at  the  center  and  3  mm  above  the  base  of  the  applicator.  The  applicator  and  the 
associated  helmet  were  placed  on  the  jig  to  simulate  the  actual  dose  delivery.  Three 
separate  measurements  were  made  for  each  applicator  to  minimize  the  effect  of  statistical 
uncertainties  from  TLD  cubes.  The  irradiation  duration  is  controlled  so  that  the  dose 
delivered  to  TLD  is  around  100  cGy  for  each  irradiation.  The  use  of  the  TLD  dosimetry 
system  follows  strictly  a  well-established  protocol  in  our  laboratory.  For  each  batch  of 
TLD,  the  sensitivity  of  each  cube,  or  chip  factor,  was  determined  by  irradiating  the  TLDs 
to  a  known  dose  and  comparing  the  TLD  readings.  To  relate  the  TLD  reading  to  dose 
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delivered  by  or  ‘®^Pd  seeds,  a  dose  calibration  was  performed  for  a  6  MV  photon 
beams  on  a  Clinac-2100C  linear  accelerator.  An  energy  correction  factor,  which  takes 
into  account  of  the  energy  response  of  the  TLD,  was  then  applied  to  yield  the  dose  given 
by  or  *°^Pd.  For  the  TLD  verification  measurements  performed  for  the  seeds,  the 
correction  factor  of  1.41  from  Meigooni  and  Nath  was  used.  Table  V  compares  the  dose 
measured  in  polystyrene  to  the  dose  calculated  to  water  at  the  same  measurement  point 
for  six  applicators.  The  measured  dose  in  polystyrene  is  on  the  order  of  5%  higher  than 
the  calculated  dose.  Relative  dose  distribution  in  selected  planes  is  being  measured  by 
using  the  GafChromic  films. 


Table  V.  Comparison  of  Dose  Measured  by  TLDs  for  Six  Applicators 


Applicator 

Measured 

(cGy) 

Calculated 

(cGy) 

Relative 
Difference  (%) 

Standard 

Deviation 

1 

124.3 

114.2 

1.089 

0.02 

L 

102.1 

96.7 

1.056 

0.00 

F 

87.9 

83.6 

1.051 

0.04 

G 

79.7 

74.8 

1.067 

0.02 

B 

87.1 

82.5 

1.050 

0.09 

H 

82.0 

78.3 

1.047 

0.04 

Experiments  were  conducted  using  seeds  with  an  initial  dose  rate  of  8  cGy/hr. 
The  response  of  the  tumors  to  treatment  was  analyzed  by  measuring  the  tumors  twice 
weekly  until  each  tumor  has  reached  a  maximum  volume  of  1  cm^  or  imtil  the  tumor  has 

regressed  and  the  animal  has 
been  free  of  tumor  for  100 
days.  Some  of  the  results  are 
shown  in  Figure  10.  Tumor 
size  growth  curve  includes 
both  the  growth  in  cell 
population  and  the  increase  in 
tumor  growth  supporting 
matrix  such  as  blood  vessels. 
Therefore  a  direct 
comparison  between 
calculation  must  include  the 
both  factors  in  the  model 
calculation.  In  the  above 
figure,  the  calculation 
includes  only  the  cell-kill  and 

cell  re-population.  The  cell 
doubling  time  in  the  tumor 
system,  as  opposed  to  the 
tumor  size  doubling  time,  was  estimated  from  the  total  tumor  cells  counted  in  the  tumor 
that  has  grown  21  days  after  an  initial  injection  of  7500  cells  to  the  rat.  The  cell  doubling 
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Figure  10. 
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time  is  estimated  to  be  2.1  days  (as  opposed  to  2.7  days  tumor  size  doubling  time,  as 
shown  in  Figure  9).  The  a,  p,  and  repair  half-time  were  determined  directly  from  the 
measurements  performed  on  the  BAl  1 12  cells.  The  calculated  curve,  with  the  radiation 
applied  at  day  17  of  tumor  growth,  show  an  initial  reduction  of  the  total  tumor  cells.  As 
time  progress  further,  the  cell  re-population  overcomes  the  radiation  induced  cell  kill  due 
to  the  low  dose  rate  and  grow  at  an  apparent  slower  growth  rate.  Modeling  the  effect  of 
tumor  matrix  grow  is  being  pursued. 

In  vivo  studies  using  an  in  vitro  assay 

In  this  experiment,  BAl  1 12  tumor  cells  were  irradiated  in  vivo  to  graded  doses  from  2  to 
20  Gy.  Following  irradiation  at  low  dose  rate,  the  tumours  were  removed  and  their  colony 

formation  ability  was  measured 
using  our  in  vitro  assay 
techniques.  Figure  1 1  shows  a 
comparison  of  the  surviving 
fraction  as  a  function  of  dose 
calculated  for  the  BAl  112  cells 
to  that  measured  in  an  in 
vivo/in  vitro  experiment  using 
CLDRI  irradiation.  The 
solid  line  represent  the 
calculated  survival  curve  using 
the  average  initial  dose  rates 
used  in  the  experiments.  The 
open  circle  represent 
calculation  using  the  actual 
initial  dose  rate  for  each 
experiment.  The  parameters  of 
a,  P,  repair  half-time,  and 
tumor  doubling  time  were 
determined  directly  from  the 
measurements  performed  on 
the  BAl  1 12  cells  (details  are 
provided  separately). 

Dose  (Gy) 

Figure  1 1 
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KEY  RESEARCH  ACCOMPLISHMENTS 


•  A  theoretical  model  based  on  incomplete  repair  during  CLDRI  has  been  developed 
for  addressing  the  questions  raised  in  the  project.  The  BED  for  implants  with  a 
mixture  of  two  radionuclides  has  been  derived  as  an  analytical  expression.  A 
manuscript  describing  this  work  has  been  submitted  for  publication  in  Intemational 
Journal  of  Radiation  Oncology  (see  Appendix  for  details). 

•  Cell  survival  curves  for  both  and  *°^Pd  were  measured  using  monolayers  of 
Chinese  hamster  cells  in  a  petri  dish  irradiated  at  low  dose  rates  using  and  ’®^Pd 
sources.  The  dose  sparing  effect  of  7  cGy/hr  relative  to  12  cGy/hr  can  be  expressed 
by  dose  modifying  factors  of  2  ±0.6  and  1.5  ±0.5  for  ^°^Pd  and  respectively.  The 
RBEs  of  **^^Pd  relative  to  were  1.2  ±  0.4  and  2.0  ±  0.5  for  7  and  12  cGy/hr, 
respectively.  In  our  system,  the  RBE  of  '*^^Pd  at  19.7  cGy/hr  relative  to  at  7.72 
cGy/hr  is  estimated  to  be  3±1.  A  manuscript  describing  this  work  is  in  preparation. 

•  Cell  survival  curves  for  *®^Pd  were  measured  using  monolayers  of  BAl  112  cells  in  a 
petri  dish  irradiated  at  low  dose  rates  using  ’°^Pd  sources.  An  orthovoltage  x-ray 
machine  was  adapted  to  produce  nearly  monoenergetic  21  keV  photons  which 
simulates  '°^Pd  photon  energies.  Using  this  x-ray  beam,  cell  survival  curves  for  the 
BAl  112  cells  in  a  petri  dish  irradiated  at  an  acute  dose  rate  were  also  measured.  We 
are  now  attempting  to  test  our  theoretical  model  for  predicting  the  CLDRI  survival 
curves  in  vitro  for  *®^Pd  sources  from  the  acute  dose  rate  exposure  data. 

•  Cell  survival  curves  for  *^^I  were  measured  using  of  BAl  112  cells  irradiated  in  vivo 
at  low  dose  rate  of  8  cGy/hr  using  the  afterloading  rat  applicator  with  *^^I  sources.  An 
orthovoltage  x-ray  machine  was  adapted  to  produce  nearly  monoenergetic  28  keV 
photons  which  simulates  ‘^^I  photon  energies.  Using  this  x-ray  beam,  cell  survival 
curves  for  the  BAl  112  cells  in  a  petri  dish  irradiated  at  an  acute  dose  rate  were  also 
measured.  We  are  now  attempting  to  test  our  theoretical  model  for  predicting  the 
CLDRI  survival  curves  in  vivo  for  ’^^I  sources  from  the  acute  dose  rate  exposure 
data. 

•  In  order  to  produce  a  consistent  dose  distribution  to  irradiate  the  tumors  transplanted 
to  different  animals  and  to  minimize  the  radiation  exposure  to  personnel  handling  the 
radioactive  seeds,  an  afterloading  seed  applicator  has  been  designed.  Twelve 
applicators  were  fabricated.  Dose  distributions  produced  by  these  applicators  were 
calculated  and  verified  by  a  series  of  dosimetry  measurements.  Several  animals  have 
been  treated  with  CLDRI  using  these  applicators  containing  '^^I  seeds  at  9  cGy/hr. 
Tumor  growth  was  significantly  slowed  by  CLDRI,  from  a  tumor  doubling  time  of 
2.7  days  in  controls  to  13  days  in  the  treated  animals.  Further  experiments  and 
analysis  are  ongoing. 

REPORTABLE  OUTCOMES 

A  manuscript  entitled  “Biologically  effective  dose  (BED)  for  interstitial  seed 

implants  containing  a  mixture  of  radio-nuclides  with  different  half  lives”  by  Zhe  Chen, 
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Ph.D.  and  Ravinder  Nath,  Ph.D.  has  been  submitted  for  publication  in  the  hitemational 
Journal  of  Radiation  Oncology.  The  manuscript  is  attached  as  an  Appendix. 

CONCLUSIONS 

We  have  made  considerable  progress  towards  the  specific  aims  of  the  project. 
Theoretical  model  for  continuous  low  dose  rate  irradiation  using  a  mixture  of 
radionuclides  has  been  developed.  Experiments  have  been  performed  using  BAl  1 12 
tumor  cells  and  Chinese  Hamster  cells  growing  in  vitro  and  BAl  1 12  cells  growing  in 
vivo  as  solid  tumors  in  WAG/rij  rats.  Radiobiology  parameters  for  these  cells  have  been 
determined  and  used  in  the  theoretical  radiobiology  model  to  improve  our  understanding 
of  the  experimental  observations.  We  have  designed  and  fabricated  applicators  for  in  vivo 
irradiations  as  well  as  developed  the  animal  care  procedures.  We  have  performed  in  vivo 
experiments  for  tumor  growth  studies  using  the  BAl  1 12  rat  model  with  seed 
applicators.  Further  experiments  with  *”^Pd  and  mixed  radionuclides  are  in  progress.  The 
in  vivo  studies  with  the  rat  model  have  been  very  difficult  to  perform  because  of  the 
nature  of  these  brachytherapy  experiments,  which  involve  long  irradiation  times.  We 
would  focus  on  completing  the  studies  using  the  rat  model  in  the  next  year.  We  expect  to 
prepare  at  least  two  manuscripts  by  the  end  of  next  grant  year. 
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Chen,  et.al.,  BED  for  permanent  seed  implants  with  a  mixture  of  radionuclides 


Abstract 

Purpose:  To  develop  a  tool  for  evaluating  interstitial  seed  implants  that  contain  a  mixture  of 
radionuclides  with  different  half-lives  and  to  examine  the  clinical  implications  of  prescribing  to 
an  isodose  surface  for  such  an  implant. 

Methods  and  Materials:  A  linear-quadratic  model  for  continuous  low  dose  rate  irradiation  was 
developed  for  permanent  implants  containing  a  mixture  of  radionuclides.  Using  a  generalized 
equation  for  the  biological  effective  dose  (BED),  the  effects  of  cell  proliferation  and  sub-lethal 
damage  repair  were  examined  systematically  for  implants  containing  a  mixture  of  radionuclides. 
The  results  were  contrasted  with  those  for  implants  using  a  single  type  of  radionuclide.  A  head  & 
neck  permanent  seed  implant  that  contained  a  mixture  of  and  *°^Pd  seeds  was  used  to 
examine  the  clinical  implications  of  the  isodose  prescription  for  such  implants. 

Results:  An  equation  of  BED  for  implants  containing  any  number  of  radionuclide  types  was 
obtained.  For  implants  containing  a  mixture  of  radionuclides  with  different  half-lives  such  as 
and  *®^Pd,  the  dose  as  well  as  its  temporal  delivery  pattern  to  a  point  is  dependent  on  the  relative 
dose  contributions  from  different  types  of  radionuclide.  It  can  vary  from  point  to  point 
throughout  the  implant  volume.  Therefore  the  quantitative  effects  of  cell  proliferation  and  sub- 
lethal  damage  repair  are  spatially  dependent  in  such  an  implant.  For  implants  containing  a 
mixture  of  and  ‘“^Pd  seeds,  the  prescription  to  an  isodose  surface  becomes  non-unique.  If  the 
prescription  dose  was  based  on  existing  clinical  experience  of  using  *^^I  seeds  alone,  mixing 
®^Pd  seeds  with  ’^^I  seeds  would  decrease  the  cell  survival  in  such  implant.  On  the  other  hand,  if 
the  prescription  dose  was  based  on  existing  clinical  experience  of  using  '°^Pd  seeds  alone, 
mixing  *^^I  seeds  with  '®^Pd  seeds  in  a  same  implant  would  create  radiobiologically  “cold”  spots 
(i.e.  an  increase  in  cell  survival  from  the  clinical  expectation)  at  locations  where  a  major  portion 
of  prescription  dose  is  contributed  by  the  *^^I  seeds.  For  fast-growing  tumors,  these  "cold"  spots 
can  become  significant. 

Conclusions:  When  cell  proliferation  and  sub-lethal  damage  repair  are  present  during  dose 
delivery,  total  dose  alone  is  no  longer  sufficient  for  a  complete  characterization  of  an  interstitial 
seed  implant.  In  order  to  avoid  radiobiological  “cold”  spots  when  radionuclides  of  different  half- 
lives  are  mixed  in  a  permanent  implant,  the  dose  prescription  should  be  based  on  the  clinical 
experience  of  using  the  longer  half-life  radionuclide.  Biologically  effective  dose  provides  a  tool 
to  start  examining  the  radiobiological  effects  of  mixing  different  type  of  radionuclides  in  the 
same  implant. 


Key  words:  interstitial  implant,  biologically  effective  dose,  iodine-125,  palladium- 103, 
brachytherapy,  radiobiology 
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L  Introduction 

Permanent  implantation  of  encapsulated  radioactive  seeds  in  tumors  has  been  used 
widely  as  a  primary  or  adjuvant  therapy  for  treating  prostate  and  head  &  neck  cancers  [2,15,26- 
27,28,29].  At  present,  seeds  that  contain  the  radionuclide  of  or  '°^Pd  are  routinely  used  for 
permanent  interstitial  implant  [23].  The  choice  of  radionuclide  for  a  given  implant  has  been 
influenced  largely  by  the  historical  development  of  radioactive  seeds  and  by  the  clinical 
experience  accumulated  with  the  use  of  each  radionuclide  [2,28].  One  rarely  mixes  seeds  of 
different  radionuclide  type  in  the  same  implant,  due  to  the  lack  of  existing  clinical  experience  as 
well  as  the  lack  of  appropriate  tools  for  evaluating  such  an  implant.  Recently,  two  permanent 
seed  implants,  each  containing  a  mixture  of  and  '®^Pd  seeds,  have  been  carried  out  in  our 
clinic  for  patients  with  head  &  neck  cancers.  The  decision  of  using  a  mixture  of  and  *”^Pd 
seeds  in  the  same  implant  was  based  primarily  on  the  clinician’s  intuition.  The  goal  of  this  paper 
is  to  develop  a  proper  tool  for  evaluating  this  type  of  implants. 

The  need  for  such  a  tool  stems  from  the  well-known  observations  in  radiobiology  that  the 
cell  survival  in  an  interstitial  implant  depends  not  only  on  the  total  dose  delivered  but  also  on  the 
temporal  pattern  of  dose  delivery  [23].  For  example,  seeds  emit  photons  with  an  average 
energy  of  28  keV  and  its  sovnce  strength  decays  exponentially  as  a  function  of  time  (with  a  half- 
life  of  60  days).  '°^Pd  seeds,  on  the  other  hand,  emit  photons  with  lower  average  energy,  about 
21  keV,  and  its  source  strength  decays  with  a  much  shorter  half-life  of  17  days.  The  variation  in 
radioactive  decay  half-life  gives  rise  to  distinct  temporal  patterns  in  the  delivery  of  a  prescribed 
dose  using  different  radionuclides.  Figure  1  illustrates  the  dose  delivered  as  a  function  of  implant 
time  for  and  ^^‘’Pd  implants.  The  time  required  to  delivery,  for  example,  80%  of  total  dose 
(dose  to  full  decay)  differs  by  a  factor  of  more  than  three  (about  140  days  for  implants  and 
about  40  days  for  *^°Pd  implants).  Such  a  drastic  difference  in  the  temproal  pattern  of  dose 
delivery  could  result  in  very  different  clinical  responses  at  a  given  total  dose  if  the  surviving  cells 
in  the  irradiated  volume  continue  to  proliferate  and  the  sub-lethally  damaged  cells  can  be 
repaired  during  the  dose  delivery  [7-9].  To  take  this  dose  delivery  difference  into  account,  the 
dose  to  full  decay  is  commonly  prescribed  at  a  value  of  145  Gy  for  implants  which  is 
considerably  higher  than  the  125  Gy  for  ‘°^Pd  implants  for  mono-therapy  of  prostate  cancer. 
When  seeds  of  different  half  lives  are  mixed  in  the  same  implant,  the  delivered  dose  as  well  as  its 
temporal  pattern  to  a  given  point  within  the  implant  are  now  depended  not  only  on  the  spatial 
locations  of  the  seeds  but  also  on  the  type  of  radionuclide  each  seed  contains.  The  temporal 
pattern  of  dose  delivery  becomes  spatially  dependent  in  such  an  implant.  A  treatment  planning 
tool  that  can  capture  the  interaction  of  the  temporal  pattern  of  dose  delivery  in  a  volume  implant 
and  the  underlying  tissue  radiobiology  is  therefore  needed  in  order  to  properly  evaluate  such  an 
implant. 

Radiation  induced  cell-inactivation,  and  the  concept  of  biologically  effective  dose  (BED), 
has  often  been  used  as  a  surrogate  of  the  biological  responses  to  ionizing  radiation  [1,5-8].  For 
acute  irradiation,  the  linear-quadratic  (LQ)  model  has  been  used  to  link  the  cell  inactivation  with 
the  total  delivered  dose.  The  LQ  model  describes  the  radiobiological  properties  of  the  irradiated 
tissue  by  two  model  parameters  (a  and  P)[7].  It  has  been  widely  used  to  examine  the  dependence 
of  the  overall  cell-inactivation  on  the  dose  fractionation  patterns  in  external  beam  radiotherapy. 
In  protracted  irradiation,  the  time  dependent  biological  processes,  e.g.  the  continuous 
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proliferation  of  the  surviving  cells  and  the  time-dependent  repair  of  sub-lethally  damaged  cells, 
may  influence  the  effectiveness  of  different  dose  delivery  patterns.  Thames  has  extended  the  LQ 
model  to  continuous  irradiation  at  a  constant  dose  rate  by  taking  into  account  the  kinetics  of  sub- 
lethal  damage  repair  during  the  dose  delivery  [27].  Later,  Dale  extended  these  concepts  further  to 
continuous  irradiation  with  exponentially  decreasing  dose  rate,  manifested  typically  in  the  low 
dose-rate  brachytherapy  using  a  single  type  of  radionuclide  [5-6].  Dale’s  formulation, 
fundamentally  equivalent  to  the  incomplete  repair  model  [27],  takes  into  account  both  the  cell 
proliferation  and  sub-lethal  damage  repair  during  a  dose  delivery  characterized  by  exponentially 
decaying  dose  rates.  These  model  introduced  the  concept  of  biologically  effective  dose  (BED)  as 
the  dose  which  would  produce  the  cell  survival  characterized  by  the  equation  5  =  for 

course  of  fractionated  radiotherapy  or  brachytherapy  [3,7].  The  resulting  formula  for  BED, 
despite  the  simplistic  nature  of  the  models  for  the  cell  proliferation  and  sub-lethal  damage  repair, 
offers  a  theoretical  tool  for  systematic  examination  of  the  interplay  between  the  temporal  pattern 
of  dose  delivery  and  the  underlying  tissue  properties.  Ling  has  used  the  model  to  assess  the 
relative  radiobiological  effectiveness  of  implants  using  radionuclide  of  different  half  lives[20] 
and  to  examine  the  biological  effects  of  dose  heterogeneity  inherent  to  interstitial  implant  [19]. 
In  this  article,  we  show  that  Dale’s  equation  for  BED  can  be  easily  generalized  for  implants 
containing  a  mixture  of  radionuclides  with  different  half  lives.  The  potential  of  using  the 
generalized  BED  formula  as  a  tool  for  evaluating  implants  containing  a  mixture  of  radionuclides 
is  examined.  The  usefulness  of  the  generalized  BED  formula  is  illustrated  by  examining  the 
clinical  implications  of  isodose  prescription  for  a  head  &  neck  implant  containing  a  mixture  of 
and  *°^Pd  seeds. 


II.  Methods  and  Materials 

1.  Dose  Rate  Calculation 


To  simplify  discussion,  an  implant  containing  a  mixture  of  radionuclides  is  defined  as 
follows.  Assume  the  implant  consists  of  N  seeds  made  of  M  different  type  of  radionuclides,  i.e. 

M 

i=\ 

where  Nj  denotes  the  number  of  seeds  containing  radionuclide  tjqie  i.  The  instantaneous  dose 
rate  to  a  point  r  at  time  t  after  the  seed  implantation  is  given  by  [22] 

M  . 

b{r,t)  =  ^D,,ir)e-^“  (la) 

/=i 

where  (f )  is  the  initial  dose  rate  due  to  seeds  containing  radionuclide  i  and  for  point-like 


source  is  approximately  equal  to 


9  A 


\ra-r 


giifii -r\)i</>a„), 


In  Eq.(l),  the  index  i  labels  the  radionuclide  type  and  the  index  /  labels  the  individual  seeds.  S,/ 
denotes  the  air-kerma  strength  for  the  seed  /  of  radionuclide  type  i.  A,,  Xi,  (^^„  )i,  and  gi(r)  denote 


the  dose-rate  constant,  the  radioactive  decay  constant,  the  anisotropy  constant,  and  the  radial 
dose  function,  respectively,  for  the  seeds  of  radionuclide  type  i.  ro  denotes  the  reference  distance 
(usually  1  cm)  at  which  the  dose-rate  constant  was  determined.  In  this  study,  the  delivered  dose. 
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isodose  distribution,  and  dose-volume  histograms  used  in  conventional  implant  dose  evaluation 
were  computed  by  using  Eq.  (1)  with  appropriate  parameters  for  and  *®^Pd  seeds  [22]. 

2.  Linear  Quadratic  Model  for  CLDRI 

The  kinetics  of  cell  proliferation  and  sub-lethal  damage  repair  for  continuous  low  dose 
rate  irradiation  (CLDRI)  were  modeled  as  follows  [5-6].  The  cell  proliferation  was  characterized 
by  a  tumor  potential  doubling  time  Tp  so  that  the  total  number  of  tumor  cells  at  any  time  t,  N(t), 
is  related  to  the  number  of  cells  present  at  a  reference  time  to,  N(to),  by 


The  fast-growing  tumors  have  shorter  Tp  and  the  slow-growing  tumors  have  longer  Tp.  The 
tumor  potential  doubling  time  is  tumor-type  dependent  and  may  vary  from  patient  to  patient  even 
among  the  same  tumor  type.  For  squamous  cell  head  &  neck  cancer,  a  “typical”  Tp  of  5  days 
have  been  quoted  in  literature  [24].  For  prostate  carcinoma,  Tp  ranging  from  10  to  60  days  have 
been  reported  in  the  literature  [10].  It  should  be  pointed  out  that  the  use  of  Eq.  (2)  assumes  that 
the  cells  in  a  target  volume  all  proliferate  at  the  same  rate  and  cell  loss  from  the  tumor  volume  is 
negligible. 


To  model  the  repair  of  sub-lethal  damage.  Dale  has  invoked  the  assumption  that 
radiation-induced  cell  inactivation  is  caused  by  the  damage  of  two  critical  targets  in  a  cell  [5-6]. 
When  a  radiation  event  damages  only  one  critical  target,  the  cell  is  considered  sub-lethally 
damaged,  which  is  repairable.  Cell  inactivation  occurs  only  when  the  other  critical  target  is 
damaged  before  the  existing  damage  is  fully  repaired.  Dale  assumed  that  sub-lethal  damage 
repairs  exponentially  with  time,  i.e.  if  a  sub-lethal  damage  was  inflicted  at  time  to,  then  the 

probability  for  it  persisting  to  time  t  is  .  Here  the  constant  p  models  the  rate  of  repair  and 

is  related  to  the  sub-lethal  damage  repair  half-time,  ,  by 


In  2 
ju 


(3) 


Average  repair  half-times  reported  for  mammalian  tissues  vary  from  0.5  to  3  hours  [3]. 


3.  Biologically  effective  dose  for  CLDRI  with  a  mixture  of  radionuclides 


The  biologically  effective  dose  has  been  defined  as  a  dose-equivalent  that  would  produce 
the  expected  cell  survival  characterized  by  an  equation  S  =  for  a  course  of  fractionated 

radiotherapy  or  brachytherapy  [3,7].  Here  a  describes  the  tumor  cell  radiosensitivity  as  given  by 
the  LQ  model  [7].  Using  the  models  of  sub-lethal  damage  repair  and  cell  proliferation  described 
above.  Dale  derived  an  expression  of  BED  for  implants  eontaining  a  single  type  of  radionuclide, 
i.e, 

0.6937;^- 

BED  =  D(T^)RE - ^ .  (4) 
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The  first  term  on  the  right-hand-side  of  the  Eq.(4)  is  a  product  of  the  total  dose  delivered  up  to 
time  Teff  and  a  factor  RE  that  characterizes  the  relative  effectiveness  of  the  radiation.  The  second 
term  on  the  right-hand-side  describes  the  decrease  in  dose  equivalent  due  to  cell  proliferation 
during  CLDRI.  In  Eq.(4),  Teff  denotes  an  effective  treatment  time  at  which  the  rate  of  tumor  cell 
proliferation  begins  to  exceed  the  rate  of  cell  inactivation  caused  by  the  instantaneous  dose 
deposited  at  Teff  (a  situation  unique  to  brachytherapy  due  to  the  exponentially  decreasing  dose 
rate  given  by  a  radionuclide.  We  will  discuss  the  determination  of  Teff  further  in  the  next 
section).  Thus,  dose  delivered  after  the  time  Teff  produces  no  net  tumor  cell  reduction.  The 
relative  effectiveness  RE  is  given  by 


i?£:  =  l  +  2(^) 


P.  r 


a  D{T^^) 


(5a) 


with 


r  = 


bl  f  1 


(1 


)■ 


fj, 


(1 


(5b) 


•  U  -XT 

where denote  the  initial  dose  rate  and  ^>(7^^  )  =  — ^  (1  -  e  '^ ) .  According  to  Eqs.  (4)  and  (5), 

X 

the  relative  effectiveness  of  a  given  radiation  is  dependent  on  the  repair  of  sub-lethal  damages. 
When  sub-lethal  damage  is  non-repairable  (i.e.  p  =  0),  y  =  D^(Teff)/2  and  RE  =  1  +  (p/a)xD(Tefr), 
which  is  the  same  as  the  equation  in  the  case  of  acute  irradiation  of  a  dose  of  D(Teff).  If,  on  the 
other  hand,  the  sub-lethal  damage  can  be  repaired  instantly  (i.e.  p  =  oo ),  then  y  =  0  and  RE 
equals  to  unit.  In  this  case,  cell  inactivation  results  from  the  simultaneous  damage  of  two  critical 
targets  in  the  cell  by  a  single  radiation  event.  In  general,  the  RE  depends  on  the  ability  of  sub- 
lethal  damage  repair  and  on  the  half-life  of  the  radionuclide.  The  presence  of  cell  proliferation, 
on  the  other  hand,  will  always  reduce  the  overall  BED  of  a  given  radiation,  as  part  of  the  dose 
has  to  be  used  to  inactivate  the  repopulated  cells.  In  this  study,  we  generalize  the  Dale’s 
expression  to  implants  containing  a  mixture  of  radionuclide  with  different  half-lives  (see 
Appendix  for  details).  For  implants  that  contains  a  mixture  of  radiouclide  with  two  different  half- 
lives,  the  y  in  Eq.(5)  is  given  by 


t>i  [  I 


/U-Ai  I  2A,  J 


^02  f  ^  n  \  ^  1 


//“/I2  [2X 


^2 


+ 


+ 


^^01^02  ,  ^  \ _ ^  Q  _  \  1 

+  ^2  ^2'^  M  \ 

^01^02  I  1 


“  A2  [  Aj  +  A2 


(1-e 


)- 


/I,  +// 


(1-c 
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(6) 
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where  and  denote  the  initial  dose  rates  at  the  point  of  calculation  produced  by  seeds 
containing  radionuclide  type  1  and  2,  respectively;  and 

+  It  is  easy  to  shown,  by  examining  Eqs.(5)  and  (6), 

that  when  radionuclide  of  different  half-lives  are  mixed  in  the  same  implant  the  resulting  BED 
does  not  equal  to  the  simple  addition  of  the  BEDs  from  each  radionuclide  type  alone. 

III.  Results 


The  use  of  BED  for  implants  containing  a  single  type  of  radionuclide  has  been  discussed 


in  several  publications  [3-6, 19-20].  In  this  section,  we  examine  the  properties  of  BED  for 
permanent  seed  implants  with  a  mixture  of  ‘^^I  and  ‘°^Pd  seeds  and  compare  them  to  implants 


using  a  single  type  of  radionuclide.  The  use  of  BED  is  illustrated  later  by  examining  the  isodose 


- c? - C?  ^  JL  •  ] 

prescription  for  a  clinical  head  &  neck  permanent  seed  implant  containing  a  mixture  of 
’®^Pd  seeds. 


I  and 


1 )  Effects  of  sub-lethal  damage  repair 

To  examine  the  effects  of  sub-lethal  damage  repair  for  permanent  seed  implants,  let  us 
consider  first  a  special  case:  a  permanent  implant  with  no  cell  proliferation  (i.e.  Tp  =  oo  in 
Eq.(6)).  In  absence  of  cell  proliferation,  the  Teff  is  determined  by  the  time  at  which  a  treatment  is 
physically  terminated  and  is  infinite  for  permanent  implants.  Setting  Teff  to  oo  in  Eq.(6)  and 
Eq.(5),  we  have 


RE  =  1  + 


■QqI  ^  -^02  I  2  -^Ol-^Ol  (^1  ^2 

P  X^  (X^  +  ju)  ^2  (^2  +  p)  (X^  +  ^2  )(A|  +  p)(X2  +  //) 


for  implants  with  a  mixture  of  radionuclides  and 


RE  =  1  +  ^-^  (8) 

a  A  +  ju 

for  implants  with  a  single  type  of  radionuclide.  Since  the  BED  is  simply  the  product  of  the  total 
dose  and  RE,  it  is  suffice  to  examine  the  effect  of  sub-lethal  damage  repair  on  RE  alone  when  the 
prescribed  dose  is  fixed. 


In  the  limiting  cases  where  the  sub-lethal  damages  are  not  repairable  (p,  =  0)  or  can  be 
repaired  instantly  (p  =  oo)  the  total  delivered  dose  alone  would  be  sufficient  to  characterize  both 
types  of  implants.  For  sub-lethal  damage  that  repairs  with  a  finite  period  of  time,  the  relative 
effectiveness  of  the  radiation  would  fall  in  between  the  values  of  the  two  limiting  cases.  The  RE 
as  a  function  of  repair  half-time  is  plotted  in  Fig.2  for  a  total  dose  of  80  Gy  (80  Gy  was  chosen  as 
it  was  the  prescription  dose  for  the  head  &  neck  case  to  be  discussed  later).  The  curves  in  the 
upper  panel  represent  implants  using  a  single  type  of  radionuclide  (  I  or  Pd)  and  the  curves 
in  the  lower  panel  are  for  implants  using  a  mixture  of  ‘^^I  or  ^°^Pd  seeds.  The  general  trend 
exhibited  in  Fig.2,  i.e.  RE  increases  with  increasing  repair  half-time,  is  intuitively  correct  as  the 
total  cell-inactivation  is  expected  to  increase  when  sub-lethal  damage  repair  becomes  less 
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efficient  (i.e.  longer  repair  half-times).  The  overall  trend  of  RE  between  ~  ^  ^1/2^  “  °° 

is  depicted  in  the  insert  of  the  upper  panel. 

Figure  2  indicates  that  the  dose  delivered  by  using  '°^Pd  seeds  alone  would  yield  more 
cell-inactivation  as  compared  to  the  same  dose  delivered  by  using  seeds  alone  when  there  is 
repair  of  sub-lethal  damage.  This  advantage  of  using  ‘°^Pd  seeds  alone  would  increase  initially  as 
repair  half-time  becomes  larger  and  would  diminish  eventually  when  sub-lethal  damage  becomes 
irreparable.  Since  the  radiobiological  effectiveness  (RBE)  arising  from  different  LET  [16], 
secondary  electrons  generated  by  and  '°^Pd  photons,  is  not  dealt  with  in  this  paper,  the 
differences  illustrated  here  are  resulted  purely  from  the  different  decay  half-lives  of  the  two 
radionuclides.  Because  the  decay  half-life  is  shorter  for  '°^Pd  seeds  as  compared  to  the  seeds, 
the  initial  dose  rate  would  be  much  greater  for  implants  using  '^’^Pd  seeds  to  deliver  the  same 
preseribed  dose.  The  trend  illustrated  in  Fig.2  is  consistent  with  the  dose-rate  effect  observed 
over  the  years  for  protracted  dose  delivery,  albeit  the  dose-rate  is  now  time  dependent  in  the 
permanent  implants. 

When  a  mixture  of  and  **^^Pd  seeds  is  used  in  the  same  implant,  the  RE  depends  on 
both  the  repair  half-life  and  the  relative  dose  eontribution  from  the  two  types  of  radionuelide  as 
shown  in  the  lower  panel  of  Fig.  2.  By  varying  the  relative  dose  contribution  fi'om  and  *®^Pd 
seeds,  the  RE  for  the  mixed-seed  implant  can  be  varied  continuously  from  that  of  using  alone 
to  that  of  using  *®^Pd  alone.  In  sueh  an  implant,  the  BED  caleulated  by  Eq.(7)  is  able  to  provide  a 
quantitative  characterization  of  the  effeet  of  spatially  dependent  differential  dose  contribution 
from  the  and  *^°Pd  seeds.  While  Fig.2  was  plotted  for  the  total  dose  of  80  Gy,  the  qualitative 
trend  shown  here  remained  the  same  for  other  delivered  doses.  The  magnitude  of  both  RE  and 
BED,  however,  would  inerease  with  the  inerease  of  total  dose  for  both  types  of  implant. 

2)  Effects  of  cell  proliferation 

Equations  (4)  and  (6)  also  capture  the  effeet  of  cell  proliferation  in  a  protracted 
permanent  seed  implant.  To  simplify  the  diseussion,  we  assume,  for  the  moment,  that  there  is  no 
sub-lethal  damage  repair.  Setting  p  to  zero  in  Eq.(6)  yields  a  BED  as  a  fimction  of  implant  time  t, 

BED(t)  =  TD{t)^  1  +  ^- 

L  ^ 

for  implants  with  or  without  a  mixture  of  radionuclides.  The  general  conclusion  remains  the 
same  for  the  case  in  which  the  sub-lethal  damage  can  be  repaired  instantly  (in  this  ease,  Eq.(9) 
would  reduce  to  BED(t)  =  TD(t)  -  0.693t/(aTp)).  In  Eq.(9),  the  TD(t)  denotes  the  total  dose 
delivered  up  to  the  implant  time  t  and  the  second  term  on  the  right-hand-side  describes  the  effect 
of  cell  proliferation.  The  presenee  of  cell  proliferation  would  always  reduce  the  magnitude  of 
BED,  sinee  part  of  the  delivered  dose  has  to  eounter  the  repopulated  cells  [5,  20].  The  amount  of 
reduction,  0.693t/(aTp),  depends  on  both  the  proliferation  rate  (0.693/(aTp))  and  the  implant 
duration  (t).  It  is  important  to  note  that  if  the  proliferation  rate  remains  the  same  during  the  dose 
delivery,  there  would  exist  a  time  at  which  the  cell  inaetivation  rate  induced  by  the  decaying 
instantaneous  dose  rate  would  beeome  less  than  the  proliferation  rate.  Beyond  this  time,  there 
would  be  no  net  eell-inactivation  advantage  fi'om  the  remaining  dose  delivered.  The  effective 
treatment  time,  Teff,  is  defined  as  the  time  at  which  the  cell  inactivation  rate  equals  to  the 
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proliferation  rate.  Equation  (9)  also  implies  that  a  minimum  total  dose  is  required  for  an  implant 
in  order  to  produce  a  net  cell-inactivation  {i.e.  for  BED  >  0)  if  cell  proliferation  is  present  at  the 
onset  of  the  implant. 

The  BED  calculated  at  Teff  depends  on  both  the  cell  proliferation  rate  and  the 
radionuclide's  decay  half-life.  Figure  3  shows  the  BED  as  a  function  of  the  potential  tumor 
doubling  time  for  a  total  delivered  dose  of  80  Gy.  Note  that  the  magnitude  of  BED  much  lager 
than  the  total  delivered  dose  due  the  assumption  that  the  sub-lethal  damage  repair  is  not 
repairable.  The  qualitative  properties  of  Fig.  3  remain  the  same  for  other  sub-lethal  damage 
repair  capabilities.  The  BED  for  implants  with  or  *°^Pd  alone  are  shown  in  the  upper  panel 
and  the  BED  for  implants  with  a  mixture  of  and  *°^Pd  seeds  are  shown  in  the  lower  panel.  As 
shown  in  Fig.3,  the  BED  increases  with  Tp  for  both  types  of  implants.  This  is  because  large  Tp 
represent  slow-growing  tumors  which  give  rises  to  small  cell  proliferation  rate,  hence  needing 
less  dose  to  combat  the  repopulated  cells.  For  a  given  tumor  proliferation  rate,  implants  that 
deliver  the  same  total  dose  within  a  shorter  period  of  time  will  result  in  a  larger  BED.  Therefore 
the  BED  for  a  ‘°^Pd  implant  will  always  be  larger  than  that  for  an  implant  for  the  same 
delivered  dose  at  any  given  tumor  potential  doubling  time  (see  upper  panel  of  Fig.3).  When  the 
two  type  of  seeds  are  mixed  in  the  same  implant  (lower  panel  of  Fig.3),  the  general  trend  of  BED 
as  a  function  of  Tp  remained  the  same  as  that  for  implants  with  a  single  radionuclide  type.  The 
magnitude  of  BED  at  a  given  Tp,  however,  becomes  dependent  on  the  relative  contributions  of 
physical  dose  from  the  two  twes  of  radionuclide.  The  resulting  BED  can  be  turned  from  that  of 
implant  alone  to  that  of  *“Pd  implant  alone  by  varying  the  relative  dose  contributions. 

3)  Permanent  seed  implant  with  a  mixture  of  radionuclides:  a  clinical  example 

When  both  cell  proliferation  and  sub-lethal  damage  repair  are  present,  Eq.(4)  and  (6)  can 
be  used  to  calculate  the  BED  for  implants  using  a  mixture  of  radionuclides.  For  the  calculated 
BED  to  be  clinically  meaningful,  one  would  have  to  know  the  patient  specific  tumor  potential 
doubling  time  and  the  sub-lethal  damage  repair  half-time.  Given  the  simplistic  nature  of  the 
biological  model  and  the  lack  of  reliable  means  of  determining  the  patient  specific  model 
parameters  at  present,  it  would  be  inappropriate  to  treat  the  calculated  BED  as  a  quantitative 
predictor  of  the  clinical  outcome  for  a  patient  implant.  However,  the  BED  can  be  indicative  of 
the  relative  merits  between  implants  that  utilize  different  physical  implant  strategies  using  a 
consistent  set  of  radiobiological  model  parameters.  It  is  on  this  premise,  we  regard  the  concept  of 
BED  a  valuable  tool  for  evaluating  the  quality  of  permanent  seed  implants.  As  shown  in  the 
previous  sections,  it  allows  a  quantitative  characterization  of  the  interplay  between  the 
complicated  temporal  pattern  of  dose  delivery  exhibited  in  implants  with  a  mixture  of 
radionuclides  and  the  underlying  radiobiological  processes.  In  the  following,  we  illustrate  the  use 
of  Eq.(4)  and  (6)  on  a  head  &  neck  implant  using  a  mixture  of  '^^I  and  '°^Pd  seeds  by  examining 
the  implications  of  prescribing  to  a  physical  isodose  line  for  such  an  implant. 

The  patient  was  a  45  years  old  male  who  presented  with  cancer  of  larynx.  A  permanent 
seed  implant  was  performed  to  the  post  pharyngeal  wall  using  ten  *^^I  seeds  and  twenty-three 
*”^Pd  seeds.  Based  on  the  reconstructed  seed  locations  from  the  seed  localization  films,  dose  to 
full  decay  was  calculated  using  Eq.(l)  and  the  isodose  lines  were  created  and  projected  to  the 
corresponding  radiographic  films.  Fig.4  shown  the  isodose  lines  (to  full  decay),  in  the  XY  plane 
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through  the  center  of  the  implant,  resulted  from  the  seeds  alone  (top  panel),  from  the  '°^Pd 
seeds  alone  (middle  panel),  and  from  both  type  of  seeds  implanted  (bottom  panel).  The  80  Gy 
isodose  line  was  chosen  by  the  physician  for  the  implant.  The  initial  source  strength  was  0.64  U 
for  seeds  and  0.88  U  for  ’°^Pd  seeds. 

The  basic  issue  regarding  the  dose  prescription  for  such  an  implant  is  as  follows.  For  a 
permanent  implant  using  only  a  single  type  of  radionuclide,  the  total  dose,  TD,  at  any  given  point 

is  related  to  the  initial  dose  rate,  , 

TD  =  \MT^i^D^  (10) 

For  this  type  of  implant,  a  prescription  to  an  isodose  line  is  equivalent  to  a  prescription  to  a 
corresponding  isodose-rate  line.  In  other  word,  the  initial  dose  rate  is  fixed  once  a  prescription  to 
total  dose  is  established.  The  temporal  dose  delivery  pattern  is  the  same  throughout  the  implant. 

When  an  implant  contains  a  mixture  of  two  radionuclides,  equation  (10)  becomes 

TD  =  1  + 1  (11) 

The  total  dose  is  now  dependent  on  the  initial  dose-rates  produced  by  the  two  types  of 
radionuclide  at  the  point  of  interest.  According  to  Eq.(l  1),  a  prescription  to  an  isodose  line  can 

now  be  fulfilled  by  many  different  combinations  of  and  .  Therefore  the  prescription  to 

an  isodose  line  is  now  not  unique  with  respect  to  the  initial  dose  rates  of  the  two  types  of 
radionuclide.  Furthermore,  the  overall  temporal  pattern  of  dose  delivery  now  varies  throughout 
the  implant.  In  absence  of  cell  proliferation  and  sub-lethal  damage  repair,  this  non-uniqueness 
would  have  caused  no  detectable  clinical  consequences  as  the  total  dose  alone  would  be 
sufficient  for  characterizing  the  implant  (see  earlier  discussion).  However,  as  was  shown  earlier, 
when  cell  proliferation  and  sub-lethal  damage  repair  are  present,  the  effective  cell-inactivation 

becomes  dependent  on  the  relative  dose  contributions,  i.e.  and  ,  of  the  two  types  of 
radionuclide  and  their  temporal  patterns. 

To  examine  the  clinical  implications  of  prescribing  to  the  an  isodose  line  for  this  head  & 
neck  implant,  let  us  examine,  for  the  moment,  a  single  spatial  point  on  the  isodose  line.  The  total 
dose  to  this  point  is  80  Gy  and  the  contributions  from  the  '^^I  seeds  and  *®^Pd  seeds  to  the  point 
can  be  determined  from  the  seed  locations  in  the  implant.  Figure  5  plotted  the  calculated  BED  as 
a  function  of  tumor  potential  doubling  time  for  all  possible  dose  contributions  by  ’^^I  and  *®^Pd 
seeds.  A  repair  half-life  of  1.5  hr  and  a  a/p  of  10  were  used  in  the  calculation.  It  is  seen  from 
Fig.  5  that  the  BED  is  largest  if  the  entire  prescribed-dose  to  the  point  is  contributed  by  *°^Pd 
seeds  (for  any  tumor  potential  doubling  times).  At  a  given  Tp,  the  BED  decreases  as  the 
proportion  of  prescribed  dose  contributed  by  the  *^^I  seeds  increases,  and  attains  a  minimum 
value  when  the  full  prescription  dose  is  given  by  ’^^I  seeds.  The  reduction  in  BED  with  the 
increased  dose  contribution  by  ‘^^I  seeds  is  most  significant  for  fast  growing  tumors  and  becomes 
less  significant  for  slower  growing  tumors. 

The  BED  calculated  with  the  actual  dose  contributions  from  *^^I  and  *®^Pd  seeds  for  head 
&  neck  implant  at  19  different  spatial  locations  on  the  80Gy  isodose  line  is  plotted  in  Fig.6  for 
several  potential  doubling  times.  The  following  observations  can  be  made  from  Figs.  5  and  6. 
First,  for  fast  growing  tumors,  those  locations  on  the  prescribed  isodose  line  with  dose 
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contributed  primarily  by  seeds  will  have  considerably  more  effective  cell  kill  than  the 
locations  with  dose  contributed  primarily  by  the  seeds,  almost  a  order  of  magnitude 
reduction  in  BED;  Secondly,  the  difference  in  the  differential  cell  kill  become  less  significant  for 
slow  growing  tumors.  Therefore  if  the  dose  prescription  was  established  from  existing  implant 
experience  that  uses  seeds  only,  mixing  ’°^Pd  seeds  with  seeds  would  always  increase  the 
effective  cell-inactivation  for  the  same  dose  prescription.  On  the  other  hand,  if  the  dose 
prescription  was  established  from  implant  experience  of  using  Pd  seeds  alone,  mixing  I 
seeds  with  *”^Pd  seeds  in  the  same  implant  would  create  radiobiologically  “cold”  (less  cell- 
inactivation  than  expected  from  the  existing  clinical  experience)  spots  at  locations  where  more 
dose  is  contributed  by  the  seeds.  These  “cold”  spots  can  become  significant  for  fast  growing 
tumors.  The  dose  prescription  for  the  clinical  case  was  based  on  experience  of  using  seeds 
alone.  Therefore  mixing  the  *°^Pd  seeds  with  seeds  in  this  implant,  at  the  least,  did  not  reduce 
the  treatment  effectiveness  from  the  priori  clinical  expectation. 

For  a  complete  evaluation  of  the  effect  of  isodose  prescription,  BED  at  locations  other 
than  those  on  the  prescription  isodose  line  should  also  be  calculated  as  the  total  dose  and  initial 
dose  rates  contributed  by  and  '°^Pd  seeds  are  spatially  dependent.  The  basic  conclusions 
drawn  above  would  remain  the  same. 

Discussions 

Due  to  the  protracted  nature  of  dose  delivery  associated  with  the  permanent  seed 
implants,  the  presence  of  sub-lethal  damage  repair  and  cell  proliferation  in  the  irradiated  volume 
requires  new  tools  to  characterize  the  interplay  of  the  physical  dose  delivery  characteristics  and 
the  underlying  biological  processes.  As  shown  in  this  work,  the  BED  formula  derived  by  Dale 
[5-6]  can  be  easily  generalized  for  permanent  implants  containing  a  mixture  of  radionuclides. 
The  general  conclusions  drawn  from  the  model  calculation,  as  shown  in  the  previous  section,  are 
physically  plausible  and  conform  to  the  common  intuition.  However,  it  should  be  pointed  out 
that  the  quantitative  numbers  of  BED  should  be  viewed  in  the  context  of  the  limitations  inherent 
in  the  models.  We  hope  the  concept  of  BED  and  the  generalized  BED  equation  for  implants  with 
a  mixture  of  radionuclides  would  help  stimulating  more  research  interests  in  exploring  a  realistic 
description  of  the  interplay  of  physical  dose  delivery  and  the  underlying  radiobiological 
processes. 

It  should  be  emphasized  that  radiobiological  modeling  is  intrinsically  organ  specific.  For 
a  given  organ,  factors  important  to  the  calculation  of  BED  include  the  physical  dose  and  its 
temporal  dose  delivery  pattern,  tumor  intrinsic  dose  response  characteristics  (a  and  P),  potential 
tumor  doubling  time,  and  sub-lethal  damage  repair  rate.  Although  not  considered  in  the  current 
model,  other  factors  that  affects  the  radiobiological  responses  of  tissues  such  as  the  presence  of 
hypoxic  cells  [12],  cell  cycle  effects  [13],  and  radiation  induced  apoptosis  [17-18]  and  organ 
architecture  should  be  considered  as  well  for  a  complete  radiobiological  characterization.  At 
present,  however,  not  all  factors  are  well  understood  and/or  quantitatively  characterized.  In 
addition,  some  of  these  factors  are  known  to  vary  from  patient  to  patient  even  in  the  same  type  of 
organ  tissue.  From  this  point  of  view,  one  should  be  cautioned  in  using  these  indices  in 
quantitatively  in  dose  prescription:  they  are  not  meant  to  and  should  not  be  used  as  an  absolute 
index  for  the  treatment  outcome  on  a  given  patient.  Nonetheless,  BED  provides  the  best  available 
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description  of  the  interplay  of  physical  dose  with  the  radiobiology.  When  a  given  set  of 
radiobiological  parameters  are  used  consistently,  they  are  very  usefiil  in  estimating  the  treatment 
efficacy  of  different  brachytherapy  applications  and  in  evaluating  competing  treatment  plans. 

Dose  distribution  from  interstitial  brachytherapy  is  inherently  inhomogeneous.  Since  the 
BED  is  determined  by  the  dose  and  its  temporal  delivery  pattern  at  each  spatial  point,  calculation 
of  BED  at  each  spatial  point  is  needed  for  a  complete  evaluation  of  a  planned  implant.  With  the 
advance  of  medical  imaging  modalities  and  image-based  treatment  planning  [21],  we  hope  that 
three  dimensional  tissue  structures  can  be  identified  more  accurately  and  a  3-D  distribution  of 
BED  can  be  calculated  with  confidence  [14].  With  the  BED  distribution,  radiobiologically 
significant  "cold"  spots,  due  either  to  “cold”  spots  in  total  dose  or  to  the  "cold"  spots  in  dose-rate, 
can  be  identified  during  the  planning  of  an  interstitial  seed  implant. 

IV.  Conclusions 

A  generalized  Dale  equation  for  biologically  effective  dose  is  presented  for  interstitial 
implants  containing  a  mixture  of  radionuclides  of  different  half-lives.  It  was  shown,  based  on  the 
BED  model,  that  the  effective  cell  survival  fi-om  an  interstitial  seed  implant  is  dependent  not  only 
on  the  total  delivered  dose  but  also  on  the  temporal  pattern  of  the  dose  delivery  when  cell 
proliferation  and  sub-lethal  damage  repair  are  present.  For  implants  containing  a  mixture  of 
radionuclides,  the  dose  to  a  point  as  well  as  its  temporal  delivery  pattern  is  determined  by  the 
relative  dose  contributions  tfom  different  types  of  radionuclide  that  can  vary  throughout  the 
implant  volume.  Therefore  the  quantitative  effects  of  cell  proliferation  and  sub-lethal  damage 
repair  are  spatially  dependent  in  such  an  implant.  Biologically  effective  dose  provides  a  tool  to 
begin  examining  the  radiobiological  effects  of  mixing  different  ^e  of  radionuclides  in  a  same 
implant.  The  model  calculation  suggests  that  adding  ‘°^Pd  to  implant  would  increase  the 
effectiveness  of  cell-inactivation  while  the  opposite  is  not  true  if  the  dose  prescription  was  based 
on  the  clinical  experience  established  with  the  single  radionuclide-type  implants.  It  is  hoped  that 
this  work  would  stimulate  further  research  interest  to  improve  the  radiobiological  modeling. 
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Appendix: 

Eq.(4)  can  easily  be  derived  following  the  work  of  Dale  [5-6]  for  implants  containing  a 
single  type  of  radionuclide.  The  reader  is  referred  to  Dale’s  original  articles  [5-6]  for  the  detailed 
steps  and  the  rationale  underlying  the  derivation.  The  main  derivation  steps  are  outlined  below. 

Dale  has  invoked  the  assumption  that  a  lethal  radiation  damage  is  caused  by  the  damage 
of  two  critical  targets  in  a  cell.  When  the  two  critical  targets  are  damaged  simultaneously  by  a 
radiation  event,  the  resulting  lethal  damage  is  termed  type-A  damage.  When  a  radiation  event 
damages  only  one  critical  target,  the  cell  is  considered  sub-lethally  damaged,  which  is  repairable. 
In  the  later  case,  a  lethal  damage  results  when  the  second  critical  target  is  damaged  before  the 
existing  sub-lethal  damage  is  fully  repaired.  This  type  of  lethal  damage  is  termed  type-B  damage. 
The  type-A  damage  is  always  proportional  to  the  total  delivered  dose  irrespective  of  dose  rate 
while  the  type-B  damage  is  dose-rate  dependent  since  the  independently  damaged  targets  may  be 
repaired  over  time.  The  overall  radiation  induced  cell  inactivation  is  therefore  proportional  to  the 
sum  of  type-A  and  type-B  damages.  The  main  goal  in  deriving  the  biologically  effective  dose  is 
to  determine  the  total  probability  of  type-A  and  type-B  damages  over  the  time  period  of  a  given 
dose  delivery. 

For  an  implant  with  a  mixture  of  radionuclides,  the  instantaneous  dose  rate  to  a  point  of 
interest  at  time  t  after  the  implant  is  given  by  Eq.(3)  in  the  main  text  as 


/=1 


(Al) 


where  the  summation  is  over  different  types  of  radionuclides  used  in  the  implant.  Using  this 
expression  for  instantaneous  dose  rate,  the  derivation  for  BED  follows  almost  exactly  the  same 
steps  outlined  in  the  Appendix  of  Dale’s  article  [5-6]. 


The  type-A  damage  is  proportional  to  the  dose  delivered.  Therefore,  the  total  type-A 
damage  accumulated  over  a  period  of  time  T  after  seed  implantation  is  given  by 


Total  type  -  A  damage  =  — —  (l-e  )  (A2) 

I 

The  type-B  damage  is  determined  by  the  joint  probability  of  two  targets  being  damaged  at 
different  times  by  two  different  radiation  events.  By  assuming  the  sub-lethal  damage  repairs 
exponentially  with  time,  the  total  type-B  damage  accumulated  over  the  same  period  of  time  T 
after  seed  implantation  can  be  determined. 


Total  Type  -  B  damage  = 


1 


//-A,. 


(l-e 


~(Ai+Aj)T .  ^ 


Aj  +// 


(l-e 


-(jU+A.j  )T , 


(A3) 


In  the  above  expressions,  a  and  P  are  the  linear  and  quadratic  coefficients,  respectively,  of  the 
linear-quadratic  cell  inactivation  model,  p  is  the  repair  time  constant  that  characterizes  the  rate  of 
sub-lethal  damage  repair. 
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With  Eqs.(A2)  and  (A3),  the  ratio  of  total  lethal  damage  to  the  type-A  lethal  damage 
alone,  termed  as  relative  effectiveness  by  Dale,  can  be  determined  as 


RE  =  1  + 


Total  type  -  B  damage 
Total  type  -  A  damage 


(A4) 


EZ  , 

a 


_ l_n . 

^  ^  ^  Xj+p 


-(M+^j)T 


X:  +  A; 


) 


Equation  (A4)  is  the  main  result  of  this  derivation.  The  biologically  effective  dose  (BED) 
for  an  implant  with  a  mixture  of  radionuclides  is  therefore  given  by 


BED  =  2; — (1  -  )RE  -  KT,j^ 

1  X; 


(A5) 


where  K  =  ln2/(aTp)  and  Tp  is  the  tumor  potential  doubling  time.  Note  that  Teff  is  the  effective 
treatment  time  at  which  the  tumor  cell-proliferation  rate  exceeds  the  tumor  cell-inactivation  rate 
caused  by  the  instantaneous  dose  rate  at  that  time. 


As  expected,  for  implants  with  a  single  type  of  radionuclide,  i  =  j  =  1,  Eq.  (A4)  and 
Eq.(A5)  reduces  to  Dale’s  original  equation 

RE  =  1  +  2(^)Z)o  .J-—!— 1±(1  -  )  -7^0 

a  n-X{\.-e  )  2X  X  +// 


-(/i+A)r 


(A6) 
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Figure  Captions; 

Figure  1  Accumulative  dose  delivered  as  a  function  of  implant  time  for  and  ^®^Pd 
implants.  Note  that  the  vertical  axis  was  plotted  as  the  ratio  of  delivered  dose  to 
the  dose  to  full  decay  for  the  two  radionuclides. 

Figure  2  Relative  effectiveness  calculated  as  a  function  of  repair  half-time  for  implants 
using  a  single  type  of  radionuclide  (upper  panel)  and  for  implants  using  a  mixture 
of  and  **^^Pd  seeds  (lower  panel).  The  insert  in  the  upper  panel  illustrates  the 
behavior  of  RE  over  a  large  range  of  repair  half-time. 

Figure  3  Biologically  effective  dose  calculated  as  a  function  of  potential  doubling  time  for 
implants  using  a  single  type  of  radionuclide  (upper  panel)  and  for  implants  using 
a  mixture  of  and  *®^Pd  seeds  (lower  panel).  Sublethal  damages  are  assumed 
not  repairable  in  this  plot.  See  text  for  a  discussion  on  the  magnitude  of  the  BED 
in  this  case. 

Figure  4  Dose  distribution  on  the  XY  plane  for  the  head  &  neck  implant  using  seeds 

alone  (top),  ’°^Pd  seeds  alone  (middle),  and  from  the  mixed-seed  implant 
(bottom).  The  80  Gy  isodose  line  on  the  bottom  panel  was  chosen  as  the  clinical 
prescription  isodose. 

Figure  5  Biologically  effective  dose  for  an  80  Gy  implant  using  a  mixture  of  I  and  Pd 

seeds,  as  a  function  of  tumor  potential  doubling  time  Tp  and  the  fraction  of  total 
dose  delivered  from  the  seeds.  A  sub-lethal  damage  repair  half-life  of  1.5  hr 
and  a  a/p  of  10  for  tumor  were  used  in  the  calculation. 

Figure  6  Biologically  effective  dose  at  19  spatial  locations  along  the  80  Gy  prescription 

isodose  line  for  the  head  and  neck  implant  with  and  **^^Pd  seeds  for  various 

tumor  potential  doubling  times.  A  sub-lethal  damage  repair  half-life  of  1.5  hr  and 
a  a/p  of  10  for  tumor  were  used  in  the  calculation. 
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